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ABSTRACT 

CCBi^normalized flat (matter plus cosmological constant) and open Cold Dark Mat- 
ter (CDM) models are tested by comparing their expected Hubble flow variations and 
the observed variations in a Type la supernova sample and a Tully Fisher cluster 
sample. The test provides a probe of the CDM power spectrum on scales of 0.02ft, 
Mpc -1 ^ k ^ Q.2h Mpc" 1 , free of the bias factor b. The results favor a low matter 
content universe, or a flat matter-dominated universe with a very low Hubble constant 
and/or a very small spectral index n ps , with the best fits having f^o ~ 0.3 to 0.4. The 
test is found to be more discriminative to the open CDM models than to the flat CDM 
models. For example, the test results are found to be compatible with those from the 
X-ray cluster abundance measurements at smaller length scales, and consistent with 
the galaxy and cluster correlation analysis of Peacock and Dodds (1994) at similar 
length scales, if our universe is flat; but the results are marginally incompatible with 
the X-ray cluster abundance measurements if our universe is open. The open CDM 
results are consistent with that of Peacock and Dodds only if the matter density of 
the universe is less than about 60% of the critical density. The shortcoming of the test 
is discussed, so are ways to minimize it. 
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1 INTRODUCTION 

In a previous paper on Hubble flow variations (Shi 1997a), 
I calculated the limit on the power spectrum shape param- 
eter r for flat CDM models. In this paper I extend the test 
of Hubble flow variations to open CDM models, and do a 
likelihood analysis of both type of models. Furthermore, I 
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will compare the results to those from other methods that 
test the CDM power spectrum. Finally, I will discuss the 
major shortcoming of this test, and possible improvements 
that can be made to minimize it. 

Before going into details, I would like to point out first 
that the method outlined here is different from that of Jaffe 
and Kaiser (1995). In Jaffe and Kaiser (1995), multi-mode 
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deviations from a pure Hubble expansion, including the bulk 
motion and shears were investigated for samples as a whole. 
In the method presented here, the variations of the Hub- 
ble expansion, corresponding to the isotropic component of 
Jaffe and Kaiser's multi-mode deviation, are investigated 
within samples for every subsample (excluding those with 
too small sizes such that non-linearity becomes an issue). 
Therefore the Hubble flow variation method is sensitive to 
the scale dependence and the shape of the variations, while 
the method of Jaffe and Kaiser is not. The two methods 
are however based on the same theoretical premise: density 
fluctuations give rise to peculiar velocities through gravity, 
and therefore from the amplitude and direction of peculiar 
velocity fields one can infer the underlying density fluctu- 
ations. It will be interesting to combine the two methods, 
by testing the variations of all modes of peculiar velocities, 
though it will be significantly more cpu intensive. 

The formalism of testing models using Hubble flow vari- 
ations has been reviewed in Shi (I997a,b). Here a summa- 
rization suffices. The deviation from a global Hubble expan- 
sion rate Ho, SH (= H^ oca i — Ho), of a sample is 



g H - n \ ^1 Uir i r l 



- B' 1 q q ~ 
^ cr. 



(1) 



where its bulk motion is 



and 

q q q' q q 

In the equations, r q is the position of object q in the sample 
(with earth at the origin), and S q (= cz q — Hor q ) is its 
estimated line-of-sight peculiar velocity with an uncertainty 
<j q . Spatial indices i, j run from 1, 2, 3, and identical indices 
indicate summation. 

Formally SH can always be expressed in the form 



S q r q r q ir 3 , 



i q' 



r q f q r q ,f J q , 



5H = d 3 rW(r)S(r) 



(4) 

where W (r) is the window function of the SH measurement, 
and S(r) is the estimated line-of-sight peculiar velocity field. 
It is not hard to see from eqs. (|l]), (0) and Q that 

W(r) = S-ij^I^tr-r,)- 

i 

(A-RB- 1 )^^^-^)- 

i 

q q ' i q " H 

If most of the uncertainty in S q comes from the uncertainty 
in measuring the distance r q , which is true at scales beyond 
~ 5000 km/sec, then o q tx r q . As a result, the window func- 
tion W scales linearly with Ho- 

The estimated line-of-sight peculiar velocity of object 
q, S q , is related to its true peculiar velocity v(r q ) by 



Sq = Vi{r q )f q + 6 q , 



(6) 



where e q is the uncertainty of the estimate with the standard 
deviation of a q . Therefore, SH can be broken into two parts: 
the true deviation 8H^"' and the noise SH^\ Eq. (§) then 
becomes 



d 3 rW i (r)v i (r)+ / d 3 r W(r) e,(7) 



SH = 5H (v) +5H (c) : 

where W l (r) — W(r)r l . The variance of 5H (v) depends on 
the density power spectrum of our universe P(k), the matter 
content of our universe Qo (the matter density divided by 
the critical density), and the global Hubble constant Ho, in 
the following way (Shi 1997a, b): 

|2P(fc) 



((SH M Y) = H o^o 2 / d 3 k\W\k)k 



k 2 



where the Fourier transform of W"(r), W(k), is 

wUk) = B 1 rv^-e' kr " 

W (2tt)3/2 [2^ of 

i 
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(9) 



q q' * q" * 

The variance of 5H^ depends on the sample measures and 
_ffo in the following form: 

((SH (C) ) 2 ) =Hl{B- 1 +B- 2 {A-RB- 1 )~ l 1 R a ). (10) 

If the value of Ho is precisely known, the above vari- 
ances can be directly calculated for a sample, and be com- 
pared with the observed deviations. But the value of Ho 
is still controversial, which leaves us the only option of in- 
vestigating the relative variation of Hubble flows within a 
sample. In other words, if the expansion rate of a sample 
with AT objects is Hn, and the expansion rate of a sub- 
sample with n (< N) objects is H n , a comparison can be 
made between the variation 5H n N / Hn = (H n — Hn)/Hn 
and its theoretical expectation without knowing the abso- 
lute value of Ho. Under the condition that Hn ~ Ho (i.e., 
((Hn - H ) 2 /H'i) 1/2 < 1), the variance of 5H nN /H N is 
5H n 



Hn 
SHi v) x 2 



+ 



H 

SHi e) * 2 



Ho 



M#Y_ fSH. 



mm) 



5H 
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Ho J V Ho J \ Ho J V Ho 
where variances of 8H n and 8Hn are calculated through 
eqs. (0) and (jll]) using their respective window functions 
Wn(k) and W n (Il). The cross correlation of the expected 
variations is 



8H„ 



(v) 



Ho 



8H 



(v) 



Ho 



= fi,V 2 J d\Re[W„(k)kW N * (tyk^ 



k 2 : 



(12) 



where Re[...] denotes the real part of the argument, and * 
denotes complex conjugation. The correlation of the noises 



mmy-imy 

The size of the samples I discuss in this paper is sufficiently 
large that their cosmic+sampling variance ((8H n / Ho) 2 ) 1 ^ 2 
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Figure 1. Hubble flow variations (with la errorbars) vs. the 
depth of subsamples of the Type la SN sample. Curves are noise- 
free lcr expectation from three models. 
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Figure 2. The same as fig. 1 but for the TF cluster sample. 

is less than 3 percent (Shi & Turner 1998). Eq. (0) is there- 
fore a very good approximation. The calculation of the vari- 
ation 8nN and its variance, on the other hand, has taken 
full account of both cosmic variance and sampling variance. 
They therefore do not depend on the sampling size of the 
(sub)sample. Of course for a small size subsample with too 
few objects, the variance of SH u n becomes large, so that a 
comparison between observations and the theoretical expec- 
tation will yield a less significant result. But the result is 
still statistically sound. 

In figures 1 and 2 I plot the Hubble flow variation 
8H n N, compared with the noise- free model expectation of 
the standard deviation ((8H n N / Hn) 2 ) 1 ^ 2 , as a function of 
the maximal depth R of subsamples (defined to include n 
most nearby objects), for a sample of 20 Type la supernovae 
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(SNe) (Riess et al. 1996) and a sample ol 36 clusters with 
Tully-Fisher (TF) distances (Willick et al. 1997), respec- 
tively. Three representative models are chosen for compar- 
ison: (1) the standard CDM model (sCDM), with Q = 1, 
h = ifo/(100km/sec/Mpc) = 0.5 and the power spectrum 
index n ps = 1; (2) the tilted CDM model (tCDM), with 
flo = 1, h = 0.5 and n ps = 0.7; (3) the vacuum energy 
A dominated flat CDM model (ACDM), with A = 0.7, 
fio = 0.3, h — 0.7, and n ps = 1. The power spectra of 
these models have the same functional form and COBE- 
normalization as in Shi (1997a). The Type la SN sample 
does not show significant detection of Hubble flow variations. 
It therefore favors models with smaller power on ~ 40ft" 1 
to 200h~ 1 Mpc scales, because otherwise its Hubble flow 
variation would be significantly larger. The TF cluster sam- 
ple shows a significant detection of Hubble flow variations 
on 45ft -1 to 60/i _1 Mpc scale. Its implication on models, 
however, is not obvious. 

To quantify the implications of figs. 1 and 2, a likeli- 
hood analysis is needed. Since different subsamples are not 
independent, we need to know the expected correlation be- 
tween the Hubble flow variations of two different subsam- 
ples, SHnM and 8H m M (m < n). This correlation is 



f dH nN \ ( dH mN \\ 
^o 2 / d 3 ifcRe| 



ki X 



WL (k) - w J N (k) 



+ 



} 



P(k) 
k 2 



(14) 



Given vector (AH)„ = SH„ N (n = n min , n min + 1, N — 1) 
measured from a real sample D, the likelihood of a cosmo- 
logical model with a set of parameters 6 is 



P(0\DI) = 



Af|£|V2 
with a normalization 



exp 



-{AH) T (T,)- 1 (AH) 



(15) 



Af ■ 



dO 



|E|iA 



• exp 



-(AJf) T (E)- 1 (AH) 



(16) 



1 test 5H n N here instead of H n — H n -i as in Shi (1997a), 
because no advantage was found using the latter quantity. 
The statistical tests of the two quantities, however, are al- 
most equivalent and give similar results. 

2 RESULTS AND DISCUSSIONS 

I apply eqs. (0) to ([ltj ) to two sets of COBE- normalized 
CDM models, the open CDM models (O < 1 and A = 0) 
and the flat ACDM models (fio + A = 1). The parameters of 
the models are chosen to be fio, h and n pa . For open (flat) 
CDM models, I choose the prior distribution of the param- 
eters to be 0.2 < f2 < 1 (0.3 < Q < 1), 0.5 < h < 0.8 
and 0.7 < n ps < 1.2, with uniformly distributed likelihood. 
The ranges of parameters are chosen conservatively to reflect 
measurements of mass densities on the galaxy cluster scale 
(Carlberg, Yee and Ellingson 1997), quasar lensing statistics 
in the case of flat CDM models (Kochanek 1996), Cosmic 
Microwave Background (Lineweaver and Barbosa 1997) and 
the Hubble constant (Freedman 1997). The final likelihood 
distribution functions after taking into account Hubble flow 
variations are then calculated according to eq. (|l5[), with 9 
being a three dimensional parameter space (flo, h, n ps ). Fig- 
ures 3-6 show on slices of constant n ps the 68% and 95% 
C.L. contours in the 3-D parameter space, for open CDM 
models and flat CDM models, using the Type la SN sam- 
ple and TF cluster sample. Only subsamples with a depth 
of more than 40/i -1 Mpc are used to ensure the validity of 
linear perturbation theory of gravity. 

Figures 3 to 6 clearly show that models with smaller 
powers are favored by both samples. Models with flo ~0.3 
to 0.4 fit data quite well regardless of the choice of h, n ps or 
the geometry of the universe, flo = 1 models, on the other 
hand, are strongly disfavored unless n ps and h approach the 
© 1997 RAS, MNRAS 000, 000-000 
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Figure 3. The 68% (the dashed line) and 95% (the solid line) 
C.L. contours in the (f2o,/i,n ps ) parameter space on different 
ips = constant slices, for open CDM models using the Type 
la SN sample. Areas to the left of the contours are allowed at 
their corresponding confidence levels. Shaded regions are allowed 
at the 2<r level by the X-ray cluster temperature function. The 
cross denotes the best fit. (a) n ps = 0.725; (b) n ps = 0.875; (b) 
n ps = 1.025; (b) n ps = 1.175. 
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Figure 4. The same as fig. 3 but using the TF cluster sample. 
Areas sandwiched between two dashed lines are allowed at 68% 
C.L. Areas to left of the solid line arc allowed at 95% C.L. 



lower end of their allowed ranges. The result diverges from a 
previous belief that peculiar velocity fields favor a high mat- 
ter density universe. The figures also show that the Hubble 
variation method is a more discriminative test to COBE- 
normalized open models. 

Also plotted in the figures is the parameter space al- 
lowed by the X-ray cluster temperature function constraint, 
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Figure 5. The same as fig. 3 but for flat ACDM models. 
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Figure 6. The same as fig. 3 but for flat ACDM models using 
the TF cluster sample. 



a s Q% = 0.52 ± 0.04 with a = 0.52 - 0.13ft for ACDM and 
a = 0.46 - O.lOQo for open CDM (Eke et al. 1996; see also 
similar results of Viana and Liddle 1996, and Pen 1996), at 
the 2<r level, if its error is taken at a face value. Since the X- 
ray cluster temperature function probes the power spectrum 
on ~ 8/i _1 Mpc scale, while Hubble flow variation calcula- 
tion probes mainly scales from ~ 30/i _1 to ~ 200/i _1 Mpc 
(as shown in figure 7), their consistency is not automatically 
guaranteed. Figures 3 to 6 show that the two methodologies 
do give consistent results for individual samples. 

When both samples are used to calculate a joint like- 
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Figure 7. The window function vs. its scale, (a) the TF cluster 
subsamplcs with n = 15 (the solid line) and n = N — 1 = 35 (the 
dashed line, multiplied by 100); (b) the Type la SN subsamplcs 
with n = 8 (the solid line) and n = N — 1 = 19 (the dashed line, 
multiplied by 10). Only one direction of k is shown. The plotted 
function looks similar in other k directions. 

lihood of models, models with smaller powers are more 
strongly favored, as seen in figure 8 and 9. However, while 
figure 9 shows that the Hubble flow result is still com- 
patible with the X-ray cluster constraints in the COBE- 
normalized ACDM models, figure 8 indicates that the two 
are marginally incompatible in the COBE-normalized open 
CDM models. Does it hint that the flat universe with a vac- 
uum energy is a more likely scenario for our universe? It 
is probably premature to say so, because the inconsistency 
is only marginal. At this stage, it may tell us more about 
the smoother distribution of the likelihood in ACDM mod- 
els relative to open CDM models, than about the favoring 
of one class of models against the other. If, however, the gap 
between the Hubble variation result and the X-ray cluster 
result widens in open CDM models in the future, we may 
need to seriously consider its cosmological implications. 

It is also interesting to compare the Hubble variation 
result to the power spectrum inferred from the galaxy and 
cluster correlation analysis of Peacock and Dodds (1994). 
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Figure 8. The joint likelihood distribution using both samples, 
for open CDM models. Legends are the same as in figure 3 to 6. 
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Figure 9. The joint likelihood distribution using both samples, 
for flat ACDM models. Legends are the same as in figure 3 to 6. 

Since the two approaches both test ~ lOO/t -1 Mpc scales, 
certain consistency is expected. Figure 10(a) plots QQ 2 P(k) 
(which the Hubble flow variation truely measures) of a num- 
ber of open CDM models along the 95% C.L. contour in fig. 
8. Since these models are marginally allowed at 95% C.L., 
the upper envelope of their Oo' 2 P(fc) curves in the tested 
range of 0.02 Mpc -1 k/h ^ 0.2 Mpc -1 represents a rea- 
sonable upper limit to Q,Q 2 P(k) of the COBE-normalized 
open CDM models. A similar limit can also be obtained for 
the COBE-normalized ACDM models. Figure 10(b) shows 
that this upper limit for the COB&normalized open CDM 
© 1997 RAS, MNRAS 000, 000-000 
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models is lower than the Qo' 2 P(k) of Peacock and Dodds 
(1994), had our universe had S7q approaching 1. This up- 
per limit is consistent with Peacock and Dodds' result only 
if the universe has Qo ^ 0.6, given the fig °' 3 dependence 
of their deduced P(k). For ACDM models, the upper limit 
on Qq 2 P(Ic) is consistent with the result of Peacock and 
Dodds even if flo — 1, although the consistency at Qo — 1 is 
marginal. The comparison once again shows that small Qo 
is favored, and that the test of Hubble flow variation is more 
discriminative to the COBE-normalized open CDM models 
than to the COBE-normalized ACDM models. 

One major shortcoming of the Hubble flow variation 
test is that the likelihood distribution is not symmetric rel- 
ative to the peak probability (figs. 3 to 6 and figs. 8 and 9) 
so that it has less power discriminating against models pre- 
dicting too small Hubble flow variations than against models 
with too large Hubble flow variations. This is due to the fact 
that the noise term in eq. (^) dominates the Hubble flow 
variation when the true density fluctuations are small. Thus 
the key to increase the testing power on models with small 
P(k) (and on models with larger P(k) to a lesser degree), is 
to reduce the noise term with better distance measurements, 
and more objects in the same sample volume. Since the ob- 
served and expected 8H n N / Hn is typically a few percent at 
a depth of 5,000 to 10,000 km/sec, the noise contribu- 
tion to SHnN I 'Hn has to be ^ 1% to ensure a significant 
detection of 8H n N /Hn and a little skewed likelihood distri- 
bution. For Type la SN samples, where distance measure- 
ment errors are typically 5% and random motions due to 
local non-linearities contribute ~ 5% of recession velocities 
at 10, 000 km/sec, the number of Type la SNe has to be 
<; 50 within ~ 10,000 km/sec to reduce the noise term to 
Si 1%. Apparently the sample used here (with the number 
around 15) is not enough (as shown by the one-sided like- 
© 1997 RAS, MNRAS 000, 000-000 
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Figure 10. Panel (a) plots Qj 2 P(fc) of open CDM models that 
lie on the 95% C.L. contours of figure 8. In (b) and (c), solid 
curves are the upper limits on £7j- 2 P(fc) from the Hubble flow 
variation test (see text) for the open CDM models (panel (b)) 
and the flat ACDM models (panel (c)). For comparison are the 
results from Peacock and Dodds (1994) for f2o = 1 (solid data 
points) and Qq = 0.6 (open data points). 

lihood distribution in figs. 3 and 5). But as the number of 
Type la SNe observed increases rapidly (for instance, there 
are 25 Type la SNe below 10,000 km/sec in an unpublished 
data set of Riess), it is hopeful that within the next sev- 
eral years a lot more precise limit can be put on the power 
spectrum from both the high end and the low end. 

For TF cluster samples, the noise contribution is dom- 
inated by distance measurement errors (~ 15%). Therefore, 
there has to be about 200 clusters within 10,000 km/sec to 
significantly boost the power of the Hubble flow variation 
method. This is not easy but still hopeful with larger and 
deeper surveys of the sky, and with distance measurements 
of clusters refined to better than about 10%. 
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